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Abstract 

GW calculations with fully self-consistent G and W - based on the iterative solution of the 
Dyson equation - provide an approach for consistently describing ground and excited states on 
the same quantum mechanical level. We show that for the systems considered here self-consistent 
GW reaches the same final Green function regardless of the initial reference state. Self-consistency 
systematically improves ionization energies and total energies of closed shell systems compared to 
GqWq based on Hartree-Fock and (semi)local density-functional theory. These improvements also 
translate to the electron density as exemplified by an improved description of dipole moments and 
permit us to assess the quality of ground state properties such as bond lengths and vibrational 
frequencies. 
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Many-bodyperturbation theory (MBPT) [l| in the GW approximation of the electronic 
self-energy jj 3|, the state-of-the-art method for the spectral properties of solids 0, S], has 
steadily gained popularity for molecules and nanosystems. In addition MBPT provides a 
prescription to extract total energies and structural properties from the GW approximation 
and therefore a consistent theoretical framework for single-particle spectra and total energies. 

Due to their numerical cost, GW calculations are mostly performed perturbatively {one- 
shot GqWq) on a set of siiigle particle orbitals and eigenvalues obtained from a preceding 
density functional theory p, Q] (DFT) or Hartree-Fock (HF) calculation. This procedure 
introduces a considerable starting point dependence 8-1^, which can be eliminated by 
iterating the Dyson equation to self-consistency. The resulting self-consistent GW (sc-GW) 
framework is a conserving approximation in the sense of Baym and Kadanoff ^Ij] (i.e. it 
satisfies momentum, energy, and particle number conservation laws). sc-GW gives total 
energies 13] free from the ambiguities of the GqWq scheme, in which the results depend 
on the chosen total energy functional [l3[]. However, as in any self-consistent theory, the 
question remains, if the self-consistent solution of the Dyson equation is unique. This issue 
is fundamentally different from the initial-state dependence of GoWo- For HF 
LDA/GGA+t/ [15j] calculations, it is well known that the self-consistency cycle can reach 
local minima instead of the global minimum. Moreover, a previous sc-GW study for the Be 
atom showed that pseudopotential calculations do not produce the same final GW Green 



]J] and 



function (and corresponding ionization potential) as all-electron calculations |l6l ]. 

In this letter, we demonstrate certain key aspects of the sc-GW approximation for 
closed shell molecules, that make sc-GW attractive as a general purpose electronic struc- 
ture method. First, the iteration of the Dyson equation produces a self-consistent Green 
function that is independent of the starting point, and determines both the ground- and 
excited-state properties (quasiparticle spectra) of a given system on the same quantum me- 
chanical level. This distinguishes sc-GW from other (partially) self-consistent GW schemes 

nn 

{171. 118I|. which do not lend themselves to total-energy calculations. Moreover the uniqueness 
of the sc-GW Green function facilitates an unprecedented and unbiased assessment of the 
GW approach, which previously was masked by the starting-point dependence of GqWo- 
Second, self-consistency improves total and quasi-particle energies compared to GqWo based 
on HF or DFT in (semi) local approximations and yields good agreement with high level 
quantum-chemical calculations and photo-emission data. Third, unlike GqWq, sc-GW yields 
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FIG. 1. (Color online) Total energy (a) and ionization energy (b) of N2 at each iteration of the 
sc-GW loop for a HF (squares) and PBE (circles) input Green function, for the aug-cc-pVQZ [l^ 
basis set. The (absolute values of the) differences arising from HF and PBE initializations vanish 
exponentially for both total energy (c) and ionization energy (d). 

an associated ground-state electron density, whose quality is e.g. reflected in the improved 
description of dipole moments. All these points taken together are essential for future devel- 
opments in electronic-structure theory such as vertex functions and beyond GW approaches. 

In the GW approximation the self-energy S is deflned in terms of the one-particle Green 
function G and the screened Coulomb interaction W as 

S(r, r', co) = tl ^G(r, r', u + u')W{r, r', uy^^ , (1) 

where is a positive infinitesimal, and spin- variables are omitted throughout this paper for 
simplicity. The screened interaction is defined as W{r,r',uj) = i;(r, r')/[l — i;(r, r')xo(r, r'a;)], 
where v = l/|r — r'| and the polarizability Xo depends explicitly on the Green function 
through the equation Xo(r, r', t) = —iG{r, r', r)G(r, r', — r). The Green function can in turn 
be expressed in terms of the self-energy through the Dyson equation 

G-' = G,'-[^-vo + Avn\ , (2) 

where Gq refers to the Green function of an independent particle system in an external 
potential Wext- accounts for changes in the Hartree potential due to density differences 
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between Go and G, and vq is the exact-exchange operator in HF or the Kohn-Sham exchange- 
correlation potential in DFT. The interdependence of Eq. [1] and Eq. [2] is the very origin of 
the self-consistent nature of the GW approach. 

It is common practice, however, to solve Eq. [1] and Eq. [2] just once, in the so-called 
Go Wo approximation. If Eq. [2] is not solved, the GqWq quasi-particle energies are given 
in first-order perturbation theory as corrections to the reference eigenvalues ({e°}) 20[ as 
e<^^ = + ReME{eT) - vom. 

In this work, Eqs. [T]|2] were solved iteratively; moreover the screened interaction W is 
also updated at each iteration taking into account the full frequency dependence of the 
polarizability Xo on the imaginary axis. In sc-GW the single particle excitation spectrum is 
given by the (integrated) spectral function: 

A{uj) = -l/n I dr Urn ImG{r,r',u) . (3) 

J r'^r 

The ground-state density n(r) also follows directly from the Green function 



n(r) = -2iG(r,r,r = 0") 



(4) 



The situation is more complicated for the total energy. As alluded to above, for a given 
Green function different prescriptions exist to compute the associated total energy such 



as the Galitskii-Migdal formula 



2l| . the Luttinger-Ward 22| or the Klein j23| functional. 



The latter two are variational in the sense that they are stationary at the self-consistent 
Green function, and therefore might provide better total energies than the Galitskii-Migdal 
formula when evaluated with non-self-consistent Green functions 2J, |25|]. However, at self- 
consistency all three approaches are equivalent. For this reason, we choose the Galitskii- 
Migdal formula 

^Tr{[u + ho]G{u)} + E,,^ , (5) 

because it is easier to implement. Equation p can be rewritten by using the equation of 
motion for the one particle Green function 26 1 



i^GM = -^Yl (r = 0-) [t,, + + 4 + T.I 



(6) 
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FIG. 2. (Color online) Difference 



ration interaction values (-Eci) 



aetween Galitskii-Migdal total energies (Eqm) and full configu- 
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-33|, with Egu evaluated from sc-Giy (diamonds), GoWo@iiF 
(squares) and GoWo@PBE (circles). PEE total energy differences with full-CI are reported as 
crosses. The calculations were performed using the aug-cc-pV5Z basis set[l9[. 

Here t is tlie kinetic-energy operator, the Hartree potential and and S'^ are the 
exchange and correlation parts of the self-energy, respectively. The trace of Eq. [5] is expressed 



as sum over basis functions in Eq. 
performed along the imaginary axis 



5] and the frequency integration can be conveniently 
26|. 

For comparison, we also computed GoWq total energies with different starting points. 
However, as indicated above, GqWq total energies are not uniquely defined, because the 
Green function and the self-energy are never on the same level. If, for example, the Dyson 
equation is not solved. Go and Sq = GqWq enter Eq. [61 If the Dyson equation is solved, the 
resulting Gi is still inconsistent with Sq. In the following we refer to the combination of 
Go and Sq in Eq. |6]as GqWq total energy and denote the corresponding starting point with 
©starting point. 

We have implemented sc-Gl^ in the all-electron electronic structure code FHI-aims 



27 



28| . Equation [T]|l] and Eq. |6] are solved in a numerical atomic orbital (NAO) basis using 
the resolution of identity technique to treat all two-particle operators efficiently 28(|. All 
calculations are performed on the imaginary frequency axis and the spectral function is 
obtained by analytic continuation to the real frequency axis 28|. The analytic continuation 
constitutes the only approximation of our implementation of the sc-Giy method. More 
information on the computation of the screened Coulomb interaction in sc-Giy is given in 



ihe Supplemental Material 
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29| . Further details of the implementation will be given elsewhere 



In Fig. [T] we demonstrate for N2 that the sc-GW Green function provides total energies 
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(a) and vertical ionization energies (b) that are independent of the starting point. Figured] 
exphcitly illustrates this for initializing the self-consistency cycle with HF and DFT in the 
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-|33|. Figure [2] reports the difference 



Perdew, Burke, and Ernzerhof (PBE) [3J] generalized gradient approximation, but other 
initializations like the local-density approximation (LDA) or the simple Hartree approxima- 
tion produce the same final sc-GW Green function (not shown). The deviation in the Green 
function exemplified by the (absolute value of the) total energy difference (Fig. [U^c)) and 
the ionization energy difference (Fig. [Hd)) converges exponentially fast with the number of 
iterations, canceling the starting point dependence. Further tests performed on a set of 30 
closed shell molecules (see Fig. H]) confirm this fact and demonstrate that sc-GW provides 
a recipe for linking different reference systems of independent electrons (or non-interacting 
Kohn-Sham particles) to a unified interacting many-body state. 

Having established that the sc-GW solution is independent of the starting point we now 
turn to an assessment of the performance of the GW approach for ionization potentials, 
electron densities and total energies. For elements in the first two rows of the periodic table 
(i.e. Z = 1 — 10) and small molecules like H2 and LiH accurate reference data from con- 
figuration interaction (CI) calculations are available 
to these CI values for basis set converged sc-GW, GoWq@PBE and GoW^o@HF calculations. 
A subset of these has previously been calculated using sc-GW pj|] and our results are in 
excellent agreement with the published results (with a maximum deviation of ^ 0.4%). In 
ine with previous calculation for the electron gas 23, |26|, |35| , atoms and small molecules 
isl . GqWq total energies (in various flavors) tend to be too negative. The self-consistent 
treatment largely (but not fully) corrects this overestimation and provides total energies in 
more satisfying agreement with full CI. The remaining overestimation provides a clear and 
unbiased quantification of the required vertex corrections in a beyond GW treatment. 

To assess other ground state properties, like the equilibrium atomic structure, would re- 
quire atomic forces (i.e. derivatives of the total energy with respect to atomic coordinates), 
which are presently not available for sc-GW. For diatomic molecules, however, structural 
properties such as vibrational frequencies, bond lengths and binding energies can be deter- 
mined directly from the potential energy curve. Other ground state properties, e.g. dipole 
moments, can be inferred directly from the electron density also for more complex systems. 
For brevity, we only present the case of CO here and refer to a future publication for a more 



detailed discussion of ground state properties in sc-GW 
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TABLE I. Equilibrium bond length d, vibrational frequency z^vib) dipole moment fi and binding 
energy E'b of the CO dimer. Units are respectively A, cm~^, Debye and eV. All calculation were 
performed with a Tier 4 basis set. 



In Table [I we report the experimental values for the bond length d, vibrational frequency 
z/vib, dipole moment fi and binding energy E]^ of CO 36|] together with the theoretical values 
obtained from several perturbative and non-perturbative approaches. DFT in the exact- 
exchange plus correlation in the random-phase approximation (EX+cRPA) is currently re- 
garded as one of the most advanced DFT functionals. Based on PBE, it is remarkably 
accurate for the bond length and vibrational frequency of CO [37]. However, like GqWq, 
EX+cRPA exhibits a considerable starting point dependence and gives no direct access to 
dipole moments. In sc-GW, the quality of the new density, obtained through Eq. lU is 
manifested in the improved dipole moment, which is in much better agreement with the 
experimental value than PBE and HF. Additional information on the quality of the sc-GW 



electron density are reported in the Supplemental Material[29|. The vibrational frequency, 
on the other hand, is overestimated and not substantially different from the perturbative 
GqWq values. Self-consistency over-corrects the overestimation of the GqWq binding energy, 
resulting in an underestimation of about 1 eV for i?b compared to experiment. Similarly, 
the sc-GW bond length is slightly too small and is close to GoVTo^HF. This assessment of 
the GW approach for ground state properties, facilitated by sc-GW, clearly indicates where 
future challenges in going beyond GW lie. 

Finally, we turn to the description of spectral properties. For the electron gas Holm and 
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FIG. 3. (Color online) Left panel: The spectral function of benzene calculated with a Tier 2 
basis set. Vertical dashed lines are located at experimental vertical ionization energies from 38l |. 
Right panel: comparison of experimental 38| and theoretical vertical ionization energies (VIEs) 



extracted from the spectral function of benzene for sc-GM^ (diamonds), GoW^o@HF (squares) and 
GoWo@PBE (circles). 



von Barth first reported a deterioration of the spectral properties 39| in sc-GW compared 



to GoM^o@LDA calculations. For the spectra of simple solids like si 



versy then arose with some authors advocating self-consistency 40|] and others dismissing it 
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icon and sodium, contro- 



4l| . Part of this controversy can be traced back to convergence difficulties in the early 



all-electron calculations 



42| . while the influence of the pseudopotential approximation in 



GW turned out to be larger than initially anticipated [43 1. 

To test the quality of the sc-GW spectral function we chose the benzene molecule as a 
benchmark, for which the sc-GW spectral function in Fig. |3]is compared to the GoWo@iiF 
and GqWo@PBE ones calculated using Eq. [31 The vertical ionization energies (VIEs) shown 
in the right panel of Fig. [3] correspond to the peak positions in the corresponding spectra. 
All the peaks reported in the left panel of Fig. 3 correspond to occupied quasi-particle states 
and the associated energy can be directly related to ionization energies from photoemission 
spectroscopy. No evidence of collective excitations was observed for the systems considered 
in the present work. sc-GW and GqW^o^HF show the same qualitative spectral features 
being characterized by sharp (5-function like quasi-particle peaks. The GqWq spectra in 
Fig. [3] depend strongly on the starting point: HF-(PBE-)based GqWq has a tendency to 
overestimate (underestimate) VIEs. The deviation between GoWo@}iF and GqWo@PBE 
is ~ 0.5 eV for the first ionization energy and can be as large as ~ 3 eV for lower lying 
quasi-particle states (Fig. [3] right panel). Furthermore, due to overscreening GoWq@PBE 
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yields a large broadening (i.e. short quasi-particle lifetimes) for states below —10 eV. At 
self-consistency, the quasi-particle energies and lifetimes are uniquely defined, the systematic 
(over) underestimation of GqWq calculations is considerably reduced and the resulting quasi- 
particle energies are in better agreements with photoemission data 



We further assessed the qua! 
proposed by Rostgaard et al. 



ity of iic-GW quasi-particle energies for the set of 30 molecules 



44j . For brevity our results are summarized in FigJJ] and we 



refer to the Supplemental Material for the actual numerical values 29l|. In Ref. 4J] sc-GW 
was based on the frozen-core approximation, whereas our work fully accounts for the im- 
portance of core electrons in the solution of Hedin's equations [43|. Core- valence coupling is 
therefore correctly included in our implementation, and we tentatively attribute the devia- 
tion of 0.1-0.5 eV in the first VIEs between our and Rostgaard et al. 's sc-GW calculations 
to the frozen core approximation made in their work. As for benzene, GoWo@iiF tends to 
overestimate VIEs, while GoW^o@PBE underestimates. sc-GW also slightly underestimates 
the VIEs, but gives an average deviation of only 2% compared to 6% in GoWo@PBE and 4% 
in GoWo@iiF. PBE and HF present two extreme starting points. In PBE the gap between 
the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital 
(LUMO) is severely underestimated, while in HF it is considerably overestimated. This in 
part explains the behavior of GqW^o^PBE and GqW^o^HF. Since the screening strength is 
inversely proportional to the HOMO-LUMO gap, GoW^o@PBE overscreens and GoWo@}iF 
underscreens. By tuning the fraction of exact exchange in the ground states, as e.g. in 
hybrid functional, the deviation between GqWq and experiment could be further reduced 
for this data set. However, this procedure is neither predictive, nor universal or transferable, 
because different systems will require a different amount of exact exchange. To really assess 
the quality of the GW approximation, self-consistency is therefore indispensable. 

In summary, we have demonstrated that sc-GW is independent of the starting point for 
closed shell molecules. Self-consistency improves the total energy and the spectral properties 
of the test sets compared to GqWq based on HF or PBE, whereas structural properties worsen 
compared to EX+cRPA. Moreover, electron densities become accessible in sc-GW, whose 
quality is demonstrated by the improved description of the dipole moment of CO. The sc- 
GW approach therefore provides a unified theory for the electronic ground- and excited-state 
properties of many-body systems. Most importantly, sc-GW gives unambiguous reference 
data that is essential for developing vertex corrections for (bio)molecules, nanostructures and 
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FIG. 4. (Color online) First vertical ionization energy (VIE) for 30 closed-shell molecules composed 



of 2 to 8 atoms. Experimental values are taken from Ref. 
a Tier 2 basis set. 



361]. All calculation were performed using 



extended systems, in particular for the challenging class of "strongly correlated" materials. 
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